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I will discuss several topics related with the computation of structure functions in the quark model in general 
O ' and its perturbative evolution. In particular, I address this topic in the Nambu-Jona-Lasinio model of hadrons, 

^^ ' where the nucleon is constructed as a soliton. I show that the handling of the regularization procedure is crucial 

in order to obtain exact scaling in the Bjorken limit and fulfillment of sum rules. I also include some problems 

concerning the general validity of quark model calculations. 
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> ! 1. Reminder on Deep Inelastic Scattering 

oo ■ Deep inelastic scattering (DIS) provides some of the most convincing evidence for the quark sub- 

m ■ 

o 






structure |1| of hadrons. The main object of study which is experimentany measured is the so-caUed 
hadronic tensor, 

W^tip,q;s) = j- I d^xe^'^< {p,s\[r^iO,JtH0)]\p,s) . (1.1) 

, where p, s are the hadron momentum and spin respectively, J^iO is a hadronic vector or axial current 

qj I and q is the momentum transfer. Introducing the Lorentz invariants Q^ = —q^ and x = Q'^/2p ■ q, 

T^ • one has on the basis of gauge and relativistic invariance the following decomposition (omitting parity 

violating contributions) 

X 

+ {p, - 9/-^) {p. - <l/-^) ^W,{x, Q^) 



W^t{p,q;s) = (-g^u + ^)Mr,Wi{x,Q 



p-q 



+ie^ux.- gi{x, Q') + 52(x, Q') s'' - ^^p"52(x, Q') . (1.2) 



q-p 



When a spin-zero hadron is considered, as e.g. the pion, the polarized structure functions gi and 52 
are ignored. Once the hadronic tensor is computed, the form factors are given by suitable projections. 
Finally the leading twist contributions to the structure functions are obtained from these form factors 
by assuming the Bjorken limit: 

Q2^oo with x = Q^/2p-q fixed. (1.3) 

For the spin independent part the structure functions fi are the linear combinations 

MNW,ix,Q^)^hix) and ^W2ix,Q^) ^ f2ix) = 2xhix) . (1.4) 
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where the last identity is a direct consequence of the spin 1/2 nature of partons. The structure function 
is related to the quark and antiquark distribution functions, qi{x) and qi{x) respectively by 

M^) = 2 51 e^[9i(a;)+gi(x)J , (1.5) 

i=u,d,s 

with ej the quark electromagnetic charges. Actually, exact scaling is only true up to logarithms, 
log((5^/A^), due to perturbative QCD radiative corrections. The most economic way of including 
them in a calculation is by means of the DGLAP equations [^ which, using renormalization group 
arguments, relate in a linear fashion the structure functions at a given reference scale, say Qq, to the 
scale of interest, Q^, 

fi{x,Q^) = U{Q\Ql)h{x,Ql). (1.6) 

Here, U{Q'^, Ql) is a linear matrix operator, fulfilling the properties U{Qi, Q2)U{Q2, Q3) = U{Qi, Qf) 
and U{Q^,Q^) = 1 In principle, the former equation means that one only needs to know structure 
functions at a given scale, since U generates structure functions at all other scales. So far U can 
only be computed in perturbation theory. Due to asymptotic freedom, this means that in practice 
U is reliable to relate high momentum scales, although nobody really knows what high momentum 
means in the present context; one is only left to practical convergence criteria between one loop and 
two loop calculations. The form of scaling violations predicted by QCD has been tested if one relates 
experimentally measured structure functions at different scales. 

Expression ( |1.1| ) is not always convenient for practical calculations. In case the considered hadron 
represents the ground state with the specific quantum numbers, the hadronic tensor will be related to 
the forward virtual Compton amplitude 
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W^t{P,Q;s) = — Imr;^(p,g;s), (1.7) 



which is given as the matrix element of a time-ordered product of the currents 

T^Up,Q;s)=ild^Ce"'<{p,s\T{j-{O4H0)}\p,s) (1.8) 

rather than their commutator. This time-ordered product has nicer properties than the commutator, 
since it can be built by functional differentiation with respect to external sources. Applying Cutkosky's 
rules subsequently yields the absorptive part, ImT^j^, which is to be identified with the hadronic tensor 



In the parton model, an analysis of the virtual Compton forward scattering amplitude leads, in the 
Bjorken limit, to the so called quark-target scattering formula 0]. This formula is formally correct, 
but it may require some regularization, or introduction of form factors. It is obvious that not every 
choice of regularization or form factor is consistent with the original Compton amplitude one started 
with. This point is often disregarded in practical calculations. 

The calculation of structure functions at a given scale remains a challenge within QCD, and one 
may recourse to quark models to compute them [Q, g. In all considered models one is forced to 
the so-called valence quark approximation Q. This hypothesis is based on the observation that the 
momentum fraction carried by the valence quarks increases for decreasing momentum scales. So, there 
is a momentum scale, say Qq, where the valence quarks carry all the momentum. In this picture, one 
starts with only valence quarks and gluons and sea quark distributions are generated as radiative 
corrections by means of DGLAP equations. Whether or not this assumption makes sense is analyzed 



below, but one should say that it is routinely applied in practical calculations. There is another 
assumption which is usually made, related to the fact that although models of hadrons deal with 
constituent quarks, perturbative QCD evolution is applicable to current quarks. In my opinion this 
point is at present not well understood, even in models which exhibit dynamical symmetry breaking, 
i.e. which allow for both kinds of quarks, and deserves further study. 

2. The NJL model in the Pauli- Villars regularization 

For practical calculations we adopt the Nambu-Jona-Lasinio (NJL) model R of quark flavor 
dynamics. This model offers a microscopic and non-perturbative description of dynamical chiral 
symmetry breaking for low energy non-perturbative hadron physics, and is rather successful from a 
phenomenological point of view. Mesons are described as quark-antiquark excitations of the vacuum 
in terms of poles of the corresponding Bethe-Salpeter amplitude. Baryons are described in the large 
Nc limit of the model as solitons. At a formal level the NJL quark distributions can be identified 
with those of the parton model. However, a major obstacle is that the bosonized NJL model contains 
quadratic divergences which have to be carefully regularized, to comply with the symmetries of the 
system. Actually, this point turns out to be crucial to ensure the fulfillment of sum rules, i.e. weighted 
integrals of structure functions in x space which yield low energy matrix elements. 

In the NJL chiral soliton model the issue of systematically regularizing the nucleon structure func- 
tions is rather subtle and has been handled in several ways. In a first approach the contributions from 



the polarized vacuum are neglected |10| , 11 1, since polarized vacuum contributions to static nucleon 



properties turn out to be numerically small, once the self-consistent soliton has been constructed. In 
other approaches the quark-target scattering formula has been regularized a posteriori |12, 13 1, in 
a treatment which only works in the chiral limit since quadratic divergences do not appear. This 
connection to static properties for setting up the regularization description does not provide a definite 
answer for structure functions not having a sum rule; i.e., structure functions whose integral cannot 
be written as a matrix element of a local operator. 

The NJL model is defined in terms of quark fields by the Lagrangian Q 

jOn3L = q{i^ - mo)q + 2Gnjl Uq^Q? + (^^^75^)4 (2.1) 

with the quark interaction described by a chirally symmetric quartic potential. The current quark 
mass, mo, parameterizes the small explicit breaking of chiral symmetry. Using functional techniques 
the quark fields can be integrated out in favor of auxiliary mesonic fields, A4 = S + iP. According 



to the chirally symmetric interaction in the Lagrangian (2J), S and P are scalar and pseudoscalar 



degrees of freedom, respectively. This results in the bosonized action |8[ 

A^ji^ = -iNcTiAlog{i^-mo-iS + ij5P)} + ^Jd''xti[MM^] . (2.2) 

Here the 'cut-off' A indicates that the quadratically and logarithmically divergent quark loop requires 
regularization. In order to compute properties of hadrons from the action (|2.2D a twofold procedure is 
in order. First, formal expressions for the symmetry currents have to be extracted. This is straightfor- 
wardly accomplished by adding external sources with suitable quantum numbers to the Dirac operator 
and taking the appropriate functional derivatives. Ignoring effects associated with the regularization, 
the currents would be as simple as q"y^{'^^)t°'q, with t°- being the appropriate flavor generator. Sec- 
ondly, hadron states are constructed from the effective action which in turn allows one to calculate 
the relevant matrix elements of the symmetry currents. For the pion this will be a Bethe-Salpeter 



wave-function which is obtained by expanding the action (2^) appropriately in the flelds S and P. 



In the case of the nucleon we will determine a soliton configuration which after collective quantization 
carries nucleon quantum numbers [|l^ ]. 

We employ the Pauli-Villars regularization scheme since it is possible to formulate the bosonized 
NJL model Q completely in Minkowski space. This is particularly appropriate when applying 
Cutkosky's rules in order to extract the hadronic tensor from the Compton amplitude. Also, let 
us remind that in ref.|jl5|, scaling for the pion structure functions was accomplished in the Pauli- 
Villars regularization, and not in the proper-time regularization. In this context, it has been shown 



1 16] that in the Pauli-Villars regularization, unlike the more customary proper-time regularization 



where cuts in the complex plane appear |17, |^, dispersion relations are fulfilled. The Pauli-Villars 
regularization has been considered before in this context both for mesons and solitons and we refer to 
refs []T9| , pC| , ^ for more details and results in this regularization scheme. We will specifically follow 
ref |I^ because in that formulation a consistent treatment solely in Minkowski space is possible. 
The regularized action of the bosonized NJL model is then given by^ 

^NJL = Ar + Ai + — f d'^x tr [52 + p2 + 2moS] (2.3) 



Nc 



An = -i^}_^CiTTlog -BB5 + At-ie , (2.4) 

Ai = -i^Trlog[-D(D5)-'-ie] . (2.5) 

The local term in eq ( p.3| ) is the reminder of the quartic quark interaction of the NJL model. After 
having shifted the meson fields by an amount proportional to the current quark masses niQ it also 
contains the explicit breaking of chiral symmetry. Furthermore we have retained the notion of real 
and imaginary parts of the action as it would come about in the Euclidean space formulation. This 
is also indicated by the Feynman boundary conditions. When disentangling these pieces, it is found 
that only the 'real part' Ar is ultraviolet divergent. It is regularized within the Pauli-Villars scheme 
according to which the conditions^ 

2 2 

Co = 1 , Ao = , ^ ci = and ^ qA^ = (2.6) 

i=0 i=0 

hold. The 'imaginary part' Ai is conditionally convergent, i.e. a principle value description must be 
imposed for the integration over the time coordinate. A priori this does not imply that it should not 
be regularized. However, in order to correctly reproduce the axial anomaly we are constrained to leave 
it unregularized. 

Essentially we have added and subtracted the (unphysical) D5 model to the bosonized NJL 
model. Under regularization the sum, log (D) + log(D5) is then treated differently from the dif- 
ference, log (D) — log (D5). In the case of the polarized nucleon structure functions ii turns out that 
this special choice of regularization nevertheless requires further specification. 

2a. Vacuum and meson sectors 

We consider two different Dirac operators in the background of scalar (S) and pseudoscalar (P) 



fields P, PI 



iB = i^-(5 + i75P) + ^ + ^75=:«DW + ^ + ^75 (2.7) 



^We denote traces of discrete indices by "tr" while "Tr" also contains the space-time integration. 
^In the case of two subtractions we need at least two cut-offs Ai and A2. In the limit Ai — > A2 = A, we have 
E. c./(A?) = /(O) - /(A^) + AV'(A'). For instance, J2^ ^^^" = (2n - 2)A2'\ 



iD. 



-i(^ - {S - i-j^P) -f + ^75 =: ^D 



(n) 



+ ii\ 



5 • 



(2. 



Here we have also introduced external vector {v^) and axial- vector (a^) fields. As noted above the 
functional derivate of the action with respect to these sources will provide the vector and axial-vector 
currents, respectively. For later use we have also defined Dirac operators, D^'^' and Dg , with these 



fields omitted. Of course, all fields appearing in eqs ( p.TD and ( |2.8D are considered to be matrix fields 
in flavor space. It is worth noting that upon continuation to Euclidean space, D5 transforms into the 



hermitian conjugate of D |18|. 

In the vacuum sector the pseudoscalar fields vanish while the variation of the action with respect 
to the scalar field S and yields the gap equation 



1 
2G 



(m — mo) = —^JNcmy^^ Ci 



i=0 



_(£k_ 



-k"^ + m^ + A? 



l€ 



(2.9) 



This equation determines the vacuum expectation value of the scalar field (S) = m which is referred 
to as the constituent quark mass. Its non-vanishing value signals the dynamical breaking of chiral 
symmetry. Next we expand the action to quadratic order in the pion field vr. This field resides in the 
non-linear representation of the meson fields on the chiral circle 



M 



mU = m exp i — vr • r 
m 



(2.10) 



This representation also defines the chiral field U. The quark-pion coupling g will be specified shortly. 
Upon Fourier transforming to vf we find 



"4njl = 5^ 



dS 



(27r)4 
with the polarization function 



niq) ■Tx{-q) 



2Ncq^n{q^ 



1 WlQ 

2G m 



+ 



o{^^), 



(2.11) 



U{q\x) 
n(g2) 



i=0 






-k^ - x(l - x)q^ + w? + K. 



le 



and 



dxll{q ,x) , 



(2.12) 



parameterizing the quark loop. The on-shell condition for the pion relates its mass to the model 
parameters 

o 1 mo 1 



m^ 



2G m 2NcU{ml) ' 
Requiring a unit residuum at the pion pole determines the quark-pion coupling 



g'^ dq^ 



qMi 



(2.13) 



(2.14) 



The axial current is obtained from the linear coupling to the axial- vector source a^. Its matrix element 
between the vacuum and the properly normalized one-pion state provides the pion decay constant /t^ 
as a function of the model parameters 



U = ANcmgIl{ml) . 



(2.15) 



The empirical values m^ = 138MeV and f^^ = 93MeV are used to determine the model parameters. 

2b. Soliton sector 

In order to describe a soliton configuration we consider static meson configurations. In that case 
it is suitable to introduce a Dirac Hamiltonian h via 

iD^ = /3{idt - h) and iD^^^ = {-idt - h)p . (2.16) 

For a given meson configuration the Hamiltonian h is diagonalized 

h^a = e„^„ (2.17) 

yielding eigen-spinors ^a and energy eigenvalues e^. In the unit baryon number sector the well-known 
hedgehog configuration minimizes the action for the meson fields. This configuration introduces the 
chiral angle 0(r) via 

h = a-p + Pmexp if ■ T^^@{r) . (2.18) 

The eigenstates |q) of this Dirac Hamiltonian are in particular characterized by their grand-spin 



quantum number |21|. The grand-spin, G is the operator sum of total spin and isospin. Since G 



commutes with the Dirac Hamiltonian ( p.l8| ) the state exp(i7rG2)|a) is also an eigenstate of (2.18) 
with energy ea and grand-spin Ga- This rotational symmetry, which actually is a grand-spin reflection, 
will later be useful to simplify matrix elements of the quark wave-functions, ^a- 

For unit baryon number configurations it turns out that one distinct level, \I'val) is strongly bound 
Q. This level is referred to as the valence quark state. The total energy functional contains three 
pieces. The first one is due to the explicit occupation of the valence quark level to ensure unit baryon 
number. The second is the contribution of the polarized vacuum. It is extracted from the action 
( |2.3| ) by considering an infinite time interval to discretize the eigenvalues of dt- The sum over these 
eigenvalues then becomes a spectral integral p2| which can be computed using Cauchy's theorem. 



Finally, there is the trivial part stemming from the local part of the action ( |2.2D . Collecting these 
pieces we have |20, ^ 

2 

Etotm = ^(l-sign(6.ai))evai-^Ec^E|\/4 + A?-\/^F7A2J 

+mlfljd\{l-cos{Q)). (2.19) 

Here we have also subtracted the vacuum energy associated with the trivial meson field configuration 
and made use of the expressions obtained for m,,r and /,r in the preceding subsection. The soliton is 
then obtained as the profile function B(r) which minimizes the total energy Etot self-consistently. 

At this point we have constructed a state which has unit baryon number but neither good quantum 
numbers for spin and flavor. Such states are generated by canonically quantizing the time-dependent 
collective coordinates A[t) which parameterize the spin-flavor orientation of the soliton. For a rigidly 
rotating soliton the Dirac operator becomes, after transforming to the flavor rotating frame ||22[|, 



zDW = A(3 (idt -T-Vt-h^A^ and iD^^^ = A [-idt + T-Vt-h) pA^ . (2.20) 
Actual computations involve an expansion with respect to the angular velocities 

A^^A=-T-n. (2.21) 

dt 2 ^ ^ 



According to the quantization rules, the angular velocities are replaced by the spin operator 



Q 



1 



a^ 



J. 



(2.22) 



The constant of proportionality is the moment of inertia a^ which is calculated as a functional of the 

"■^ is of the order l/A'^c'.^ Hence an expansion 



soliton |22|. For the present purpose we remark that a' 

in is equivalent to one in I/Nq- The nucleon wave-function becomes a (Wigner D) function of the 



collective coordinates. A useful relation in computing matrix elements of nucleon states is [14| 



{N\U,[a^t,At,)\N) 



\mhJj\N). 



(2.23) 



3. Bjorken limit and scaling 

An important consequence of the Pauli-Villars regularization is given by the fact that in the 



Bjorken limit the model scales. A detailed analysis |23| shows that in this limit the contribution to 
the Compton scattering amplitude of the regularized Dirac determinant is given by a second functional 
differentiation of 



i=0 '^ 



"i^dW 



dW( 



.-1 



-.Q' 



(3.1) 



with respect to the vector sources. For the imaginary part of the action the expression analogous to 
(lO) reads 



A 



A,I 



.^Tr{(-DWDW)"^[Q^ 



,-1 



^d^"^+dW( 



v-l 



hQ' 



(3.2) 



In both cases, it is understood that the large photon momentum runs only through the operators in 
square brackets. With Sf^p^a = g^pgua + gpuQiicr — g^iugpa, the subscript '5' means 



Iplplv = Sp.pval'^ - iCfj.puaj''^^ whilc {jfj.Jp'yu)5 = Sppuaj'^ + ICppua^'^J^ 



(3.3) 



The fact that the sum rules enforce this extension of the regularization scheme is not all surprising 
since the derivative operator i^ fixes the Noether currents. Rather it is imposed and a consequence 
of the 'sum rules' of the model defined by D5. The D5 model, which is not physical, has been 
introduced as a device to allow for a regularization which maintains the anomaly structure of the 
underlying theory. Hence, further specification of this regularization prescription is demanded in 
order to formulate a fully consistent model. It should be stressed that this issue is not specific to the 
Pauli-Villars scheme but rather all schemes which do regularize the sum, log (D) + log (D5) but not 
the difference, log (D) — log (D5) will require the specification { \i.3\) . Since only the polarized, i.e. spin 
dependent, structure functions are effected, this issue has not shown up when computing the pion 
structure functions in the Pauli-Villars scheme. In the unregularized case (Aj = 0) the contributions 
associated to the expansion of D5 cancel in the sum (|3.lD and (|3.2|) leaving 



AM = ,^Tr 



dW 



Q' 



\-i 



(3.4) 



^When formally considering n ucleo n structure functions it will turn out that it is mostly sufficient to refer to the 
Dirac operators as defined in eqs (2.2C) rather than to explicitly carry out the expansion in the angular velocities. 



Finally, it should be mentioned that the property of scaling is not an automatic consequence of any 
regularization. For instance, in the Proper-Time regularization, there appear logarithmic corrections 
since the imaginary part of the one loop diagrams is not regularized. 

4. The regularized structure functions 



4a. Pion structure functions 

Working through the expressions, one obtains after some calculation an explicit expression for the 
hadronic tensor. For the pion we have 



d 



Mx) = Qmc9'£^[q^Uiq\x 



(4.1) 



which is trivially normalized. In the chiral limit m.,^ ^ we get 

h{x) = {^)e{x)9{l-x) 



(4.2) 



This result is surprisingly simple for it does not seem to depend strongly on any dynamical assumptions, 
and it is equivalent to set the parton distributions, Ut^{x) = d-,^{x) = 1 for vr^. In addition, if DGLAP 
evolution at LO [|l5[ and NLO [g5[ is undertaken a very satisfactory description of valence quark 
distributions in the pion as extracted in ref.|2(:] is obtained. One should say however, that the gluon 
and sea quark distributions are too steep in the low x region and too flat in the high x region |25]. In 



ref.[15|, parton distributions for the kaon are also computed. 

4h. Nucleon structure functions 

Here we will only discuss the contribution of the polarized vacuum to the nucleon structure func- 
tions on a formal and conceptual level, since up to now the final formulas have not been numerically 
evaluated. The contribution of the distinct valence level, which is not effected by the regularization, 
cf. eq ( ^.19|) , has previously been detailed Q, |l^ and numerically evaluated. After some calculation, 
we obtain for the nucleon l2J 



W^,{q) 



-.^/^E/^Vi^*"'"-^ 



(4.3) 
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[^aimhf.i7u'^M + Xes)e-'^^ 



-^a(e)QA7.^7^^a(e-Ae3)e^-^- /+(, 



[UJ, 



pole 
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-^aim\{l.i7^.h^a{i+Xes)e'^'']f-{co) 



pole 
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with the light-cone vector n = (1,0,0, 1) and the spectral functions 
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/a^H=E^^ 



UJ±€r. 



iO±en 



i=0 



Ijj^ 



el - A? + ie w^ 



ei + ie 



(4.4) 



and "pole" means computing the residue contribution of the poles of the spectral functions. The 
previous equations are similar to the decomposition into quark and anti-quark distributions. In the 
former expression Eq.([4.4|), the sum is over the continuum Dirac spectrum. In the unregularized case, 
fa{io) = while /^('^)|poie — — 47rz(^(c(; — eQ,). Apparently the hadronic tensor then indeed becomes 
a sum of quark and anti-quark distributions. However, in the Pauli-Villars regularized scheme, we 
have additional contributions from quark and anti-quark distributions with dispersion relations which 
also contain the cut-offs, Aj. Hence they differ from those dispersion relations naively expected from 
the solutions of the Dirac equation ( ^.17 ). Prom this expression structure functions can be directly 



obtained by contracting with appropriate projectors. 

Note that within the Bjorken limit the Callan-Gross relation, f2ix) = 2xfi{x), is automatically 
fulfilled. The unpolarized structure function fi{x) then becomes 



Nc fduj^^f ,3^ f dX iM.rxX f^^ ^ + ^c 



/■w - --«^/^E/.^VI^''""i^ 



«=0 " * / pole 

a \j=0 " « / pole 

X / d'( {*J,(0 (1 - as) ^a(C"'+Ae3)e-^"^ - ^UO (1 " "s) ^a(e- Ae3)e'"^} 
to leading order in I/Nq- The structure function which enters the Gottfried sum rule [pO|] should 



not be regularized in contrast to previous studies. In ref IQ this structure function has been treated 



analogously to the one of neutrino nucleon scattering associated with the Adler sum rule. As discussed 



in ref. 1 23 1 the latter indeed undergoes regularization. This example clearly exhibits that obtaining the 
formal expressions for structure functions from the defining action is unavoidable in cases when there 
is no relation to a static nucleon property. Similar expressions to leading order in 1/Nc for the 
polarized structure functions, gi{x) and g2{x) can also be found in ref. [^]. Since the soliton breaks 
translational invariance, these structure functions extend out of the interval < x < 1. The issue of 



restoring the proper support has been addressed in ref.|31| by going to the infinite momentum frame. 
Sum rules relate moments of the structure functions to static properties of hadrons and a con- 
sistently formulated model is required to satisfy these sum rules. Static properties are obtained by 
computing matrix elements of the symmetry currents. In ref.p3|, it has been shown that the prescrip- 



tion (3^) complies with the required sum rules for polarized and unpolarized structure functions. 



5. Phase-Space considerations 

The result, that the pion distribution function becomes unity in the chiral limit, seems to be a 
rather remarkable one, because it looks independent of any detailed dynamics. In addition, if a DGLAP 



evolution for the valence part is undertaken to leading order |15| and NLO [Eq] the agreement with 
experimental parameterizations is striking. The interesting aspect here is that within our calculation 
we are treating the pion as a composite qq Goldstone boson. We will derive the same result by a 
rather crude method. The probability F{x) (normalized to unity) that a quark in a hadron with N 
constituents has a momentum fraction between x and x + dx is given by 

N ^ N 

F{x) = dxi... dxN^C^ Xi - l)f{xi, ..., xn)j;t Y1 ^(^ ~ ^^) (^-1) 

1=1 i=l 



GRV-LO 

GRV-NLO 

LO 

NLO 




Figure 5.1: Valence quark distribution function in the Nucleon xV{x) = x[u{x) + d{x) — u{x) — d{x)]/3 
computed at LO and NLO at Q"^ = 4GeV^. The labels LO and NLO correspond to a initial distribution 
function xV{x) = 2x(l — x) obtained by simple phase space considerations and DGLAP evolved to the 
scale Q^ = 4GeV^ (see main text). The labels GRV-LO and GRV-NLO are the "experimental" analysis of 
ref. p4ll at the same scale. 



where f{xi, 
choice f{xi, 



, xj\f) is the probabihty distribution of the N quarks in the hadron. If we make the 
, xn) = 1, then we find 



F{x) = {N-l){l-x) 



N~2 



(5.2) 



We see that in the case A^ = 2, we have Ft^(x) = 1, which was the result we obtained previously 
for the pion. It is of course very tempting to look at the nucleon case {N = 3), for which we get 
Fn{x) = 2(1 — x). After DGLAP evolution we get a reasonably good description of the valence part 
V{x) = [u{x) + d{x) — u{x) — d{x)\/2> (see figure), much better than what it is obtained in much 
more sophisticated calculations ^]. The agreement is not surprising if one looks at the low energy 
parameterizations of GRV 95' in ref. |2^], where one sees that as x ^ 1 the valence contribution 
behaves as ~ (1 — x). Actually, after evolution one finds that for Q^ = 4GeV^ the distribution 



function behaves as F{x) ~ (1 — x) '^ which only agrees very roughly with the counting rules 
found long ago [^], F{x) ~ (1 — x)^*^"^, with M = 2 for the nucleon. It is fair to say that the 
scale at which these counting rules were derived was never made clear, and that structure functions 
are not experimentally known for x > 0.75. It is also fair to say that the analysis above was already 
envisaged in the original paper of Bjorken and Paschos [^] although to that time perturbative QCD 
evolution equations were not known, so they tried to build the full Q^ dependent structure function 
by making a judicious choice of /(xi, . . . ,xiy). Although, the results found in this section should be 
analyzed with care, it seems that the bulk of the valence distributions is understood in terms of the 
number of constituents only. Nevertheless, in both cases the gluon and sea distributions are too steep 
at low X and too flat at high x [p^ j. In addition, the /2 is reasonably well described at LO and NLO 
at Q^ 



5GeV'2 between 0.3 < x < 0.75 IJ. 
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6. Can low energy models be matched to perturbative QCD 



9 



To finish the discussion let us analyze whether or not low energy models can be matched to 
perturbative QCD in a reasonable way. The traditional recipe for this kind of calculations is that 
after the low energy model structure functions are computed, a subsequent DGLAP evolution is 
undertaken. Since these models are thought not to provide a gluon or sea distribution, the initial 
condition consists of taking only the valence part, the radiative corrections due to gluons and sea 
quarks, are automatically generated by the solution of DGLAP equations. I have looked at this 
problem from a different point of view |2^. Namely, I have taken the parton distribution functions 



at very high energies, Q^ = rn^, and I have evolved them downwards in energy, until some structure 
function becomes negative, in the interval 0.01 < x < 1. This seems a natural place to stop evolution, 
since going below would violate unitarity. The relevant question to ask is whether or not around 
this point the gluon and the sea distributions are neglegeable. There is a subtlety in the procedure, 
since NLO DGLAP equations are often used in a form which is not manifestly "reversible" , i.e. the 
property U{Q\,Q2)U{Q2,Q1) = U{Q\,Q1) is not fulfilled. After properly handling this point I have 
found that the unitarity limit takes place at the scale Qq = 580MeV. At this scale one gets 

(x),al=0.45, (x)giuons = 0.36 , (x)sea = 0.19, 05 = 0.48 (6.1) 

From the point of view of perturbation theory this scale is certainly more reliable than that required 
by the valence quark model point (defined as (x)vai = 1, (ic)giuons + {x)sea. = and which requires 
as = 1.80 ), and it is still having a non-negligible gluon and antiquark distribution. Thus, there seems 
to be a gap between the low energy nucleon model and the perturbative QCD result. In perturbative 
QCD there exists a systematic expansion in powers of as- In the quark model such a systematic 
expansion seems to be missing, and therefore "improvements" of it are a bit arbitrary. What is 
needed is a better understanding on how to systematically improve the low energy model including 
gluons and antiquarks to be really able to work on the low energy side and hopefully fill in this gap. 
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